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Abstract A 2.75-approxiniation algorithm is proposed for the unconstrained 
traveling tournament problem, which is a variant of the traveling tournament 
problem. For the unconstrained traveling tournament problem, this is the first 
proposal of an approximation algorithm with a constant approximation ratio. 
In addition, the proposed algorithm yields a solution that meets both the no- 
repeater and mirrored constraints. Computational experiments show that the 
algorithm generates solutions of good quality. 
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1 Introduction 

In the field of tournament timetabling, the traveling tournament problem 
(TTP) is a well-known benchmark problem established by Easton, Nemhauser, 
and Trick [l] . The present paper considers the unconstrained traveling tourna- 
ment problem (UTTP), which is a variant of the TTP. In the following, some 
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terminology and the TTP are introduced. The UTTP is then defined at the 
end of this section. 

Given a set T of n teams, where n > 4 and is even, a game is specified 
by an ordered pair of teams. Each team in T has its home venue. A double 
round-robin tournament is a set of games in which every team plays every 
other team once at its home venue and once as an away game (i.e., a game 
held at the home venue of the opponent). Consequently, 2(n — 1) slots are 
necessary to complete a double round-robin tournament. 

Each team stays at its home venue before a tournament and then travels 
to play its games at the chosen venues. After a tournament, each team returns 
to its home venue if the last game is played as an away game. When a team 
plays two consecutive away games, the team goes directly from the venue of 
the first opponent to the venue of another opponent without returning to its 
home venue. 

For any pair of teams i, j e T, dij > denotes the distance between the 
home venues of i and j. Throughout the present paper, we assume that triangle 
inequality {dij + djk > dik), symmetry (dij = dji), and da = hold for any 
teams i, j, k ^ T. 

Denote the distance matrix (dij) by D. Given an integer parameter u >2, 
the traveling tournament problem [3| is defined as follows. 

Traveling Tournament Problem (TTP(u)) 

Input: A set of teams T and a distance matrix D = (dij). 

Output: A double round-robin schedule of n teams such that 

CI. No team plays more than u consecutive away games, 

C2. No team plays more than u consecutive home games, 

C3. Game i at j immediately followed by game j at i is prohibited, 

C4. The total distance traveled by the teams is minimized. 

Constraints CI and C2 are referred to as the atmost constraints, and Con- 
straint C3 is referred to as the no-repeater constraint. 

Various studies on the TTP have been conducted in recent years (see [H 
[TOlfT^ for detail), and most of these studies considered TTP(3) [13], which 
was recently proved to be NP-hard by Thielen and Westphal [T^. Almost all 
of the best upper bounds of TTP instances are obtained using metaheuristic 
algorithms. On the other hand, little research on approximation algorithms 
has been conducted for the TTP. Miyashiro, Matsui, and Imahori [9] pro- 
posed a (2 -I- 0(l/n))-approximation algorithm for TTP(3). Yamaguchi, Ima- 
hori, Miyashiro, and Matsui |16] proposed an approximation algorithm for 
TTP(u), where 3 < it <C Westphal and Noparlik [TS] proposec(3 a 5.875- 
approximation algorithm for TTP(m), where 3 < u. For TTP (3), the approx- 
imation ratio of [16| is the best among them. In addition, Thielen and West- 
phal [11] proposed a (1.5 -I- 0(l/n))-approximation algorithm for TTP(2). 



^ Westphal and Noparlik's paper 1151 and the conference version of the present paper [7] 
appeared in the same conference (PATAT, 2010). 



A 2. 75- Approximation Algorithm for the Unconstrained TTP 



3 



The TTP is a simplification of an actual sports scheduling problem. Some 
further simplified variants of the TTP have been studied [13]. The circular 
distance TTP and the constant distance TTP are the problems which have 
specific distance matrices. For the constant distance TTP, Fujiwara, Imahori, 
Matsui, and Miyashiro [S] proposed approximation algorithms. 

The unconstrained traveling tournament problem (UTTP) is another vari- 
ant of the TTP, in which Constraints CI through C3 are eliminated. In other 
words, the UTTP is equivalent to TTP(n — 1) without the no-repeater con- 
straint. On some actual sports scheduling problems, the atmost constraints 
(w = 3 in particular) and the no-repeater constraint are considered. However, 
these constraints are not necessarily imposed, and the UTTP is a suitable 
simplified model for some practical scheduling problems. 

Bhattacharyya [T] recently showed NP-hardness of the UTTP. Although 
the UTTP is simpler than the TTP, no approximation algorithm has yet been 
proposed for the UTTP. The method proposed in [T^ cannot be applied to the 
UTTP because the condition w <C n is necessary. The method in [9] , proposed 
for TTP(3), can be applied to the UTTP with a few modifications. However, 
this leads to a ((2/3)n -f 0(l))-approximation algorithm for the UTTP, which 
is not a constant approximation ratio with regard to n. 

In the present paper, we propose a 2.75-approximation algorithm for the 
UTTP. In addition, the solution obtained by the algorithm meets both the no- 
repeater and mirrored constraints, which are sometimes required in practice. 
This property indicates that our algorithm also works for TTP(n — 1), which 
eliminates the atmost constraints but considers the no-repeater constraint. 



2 Algorithm 

In this section, we propose an approximation algorithm for the UTTP. A 
key concept of the algorithm is the use of the circle method and a shortest 
Hamilton cycle. The classical schedule obtained by the circle method satisfies 
the property that, for all teams but one, the orders of opponents are very 
similar to a mutual cyclic order. Roughly speaking, the proposed algorithm 
constructs a short Hamilton cycle passing all venues, and finds a permutation 
of teams such that the above cyclic order corresponds to the Hamilton cycle. 

Let G = {V, E) be a complete undirected graph with the vertex set V and 
edge set E, where \V\ = n. We assume that there exists a bijection between 
the vertex set V and the set of teams T. We put the length of edge {v, u'} G i?, 
denoted by d^^i , to the distance between the home venues of the corresponding 
teams t, G T. First, we assign aliases 0,l,...,n — 1 to teams in T as follows. 

1. For each v <^V, compute Y.v'evMv} 

2. Let V* be a vertex that attains miui.gv di>«', and designate the 
team corresponding to v* as team n — \. 

3. Using Christofides' 1.5-approximation algorithm [5] for the traveling sales- 
man problem with triangle inequality and symmetry, construct a Hamilton 
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cycle on the complete graph induced hy V \ {v*}. For the obtained cycle 
(uo, ui, . . . , u„_2), denote the corresponding teams by (0, 1, . . . , n — 2). 

In the rest of this paper, we define that the set of teams T = {0, 1, 2, . . . , n — 1} 
and the vertex set V — {uq, ui, . . . , Wn-2, '^*}- We identify the vertex v„-i 
with vq (not V*) and the vertex with w„_2 (not v*). 

Next, we construct a single round-robin schedule. In the following, a "sched- 
ule without HA-assignment" refers to a "round-robin schedule without the 
concepts of home game, away game, and venue." Denote the set oi n — 1 slots 
by 5 = {0, 1, . . . , n~2}. A single round-robin schedule without HA-assignment 
is a matrix K of which (t, s) G TxS element, say K{t, s), denotes the opponent 
of team t in slot s. Let K* be a matrix defined by 



Lemma 1 |16j The matrix K* is a single round-robin schedule without HA- 
assignment. In addition, K* is essentially equivalent to the classical schedule 
obtained by the circle method. 

Then, by the mirroring procedure, we form K* into a double round-robin 
schedule without HA-assignment. More precisely, construct a matrix {K*\K*) 
whose rows are index by teams and columns are index by a sequence of slots 
(0, 1, . . . , n — 2, 71 — 1, 71, . . . , 2n — 3). So as to complete a double round-robin 
schedule, "home" and "away" are assigned to games of {K*\K*) as follows: 

— for team t G {0, 1, . . . , 7i/2 — 1}, let the games in slots 2t, 2t4-l, . . . , 7i+2< — 2 
be home games, and let the other games be away games. 

— for team t £ {n/2, n/2 -I- 1, . . . , ti — 2}, let the games in slots 2t — ii + 2, 2t — 
71 -I- 3, . . . , 2< be away games, and let the other games be home games. 

— for team ti — 1, let the games in slots 0,l,...,7i — 2 be away games, and 
let the other games be home games. 

The obtained double round-robin schedule is denoted by K^^^. Figure[T]shows 
the schedule K^^^ of 10 teams. 

Lemma 2 The double round-robin schedule K^^^ is feasible. 
Proof. {K*\K*) is a consistent double round- robin schedule without HA- 
assignment, which satisfies the mirrored constraint. We check the feasibility 
of HA-assignment to games. Teams i and j (i < j < ?i — 1) have a game at 
slot i+ j. By the rule to assign home and away to games, team i plays a home 
game and team j plays an away game at slot i-\-j. Teams i and j {i < j = n — 1) 
have a game at slot 2i, and the rule assigns consistent home/away to the teams. 
Another game between teams i and j is held at the opposite venue. □ 




s — t (mod 71 — 1) (t ^ n — 1 and s — t ^ t (mod n — 1)), 
71 — 1 (t 7^ 71 — 1 and s — t ~ t (mod n — 1)), 

s/2 (t = n — 1 and s is even), 

(s -t- 71 — l)/2 {t ~ n — 1 and s is odd). 
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Each number corresponds to the opponent and away (home) game is denoted by A (H). 



Fig. 1 The schedule K^^^ with 10 teams. 



In addition, for each m £ {0, 1, . . . , 2n — 3}, we construct a double round-robin 
schedule by rotating slots of K^^^ through m cyclically. It means that games 
of K^Yini'iTT') ("^ G {O7 1, ■ . ■ , 2n — 3}) at slot s are equal to games of -R^drr 
at slot s + m (mod 2n — 2). Obviously, all of the schedules K^^^{m) [m e 
{0, 1, . . . , 2ri, — 3}) meet both the no-repeater and mirrored constraints. Finally, 
output a best solution among ^drr(™) (w G {0, 1, . . . , 2ri — 3}). 

Here, we estimate the time complexity of the algorithm. Christofides' al- 
gorithm requires 0{n^) time to construct a Hamilton cycle on the complete 
graph induced by F \ {v*}. For the constructed Hamilton cycle, there are 
2(n — 1) possibilities to assign teams. For each assignment of teams, we con- 
sider 2n — 2 possibilities of to G {0, 1, . . . , 2n — 3}. Each double round- robin 
schedule can be evaluated in 0{n) time on average. Thus, the time complexity 
of the algorithm is bounded by O(n^). 

In the next section, we prove that the proposed algorithm guarantees an 
approximation ratio 2.75. 

3 Approximation Ratio 

In this section, we describe the proof of the approximation ratio of the proposed 
algorithm. Designate the length of a shortest Hamilton cycle on G as r. 

Lemma 3 The following propositions hold for G. 

(1) The length of any edge is hounded by t/2. 

(2) The length of any Hamilton cycle on G is bounded by nr/2. 

(3) E ^ 

vev v'ev\{v} 

(4) ^ di,y, < nr/A. 
vev\{v} 

Proof. (1) For the edges {i,j} and {j, i}, the sum of their lengths is at most 
the length of a shortest Hamilton cycle. Thus, the length of the edge {i,j} is 
bounded by r/2 with symmetry. 
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(2) This immediately follows from Property (1). 

(3) Given a shortest Hamilton cycle H = {uq, ui, . . . , w„-i) on G, let 

dui,u,+t + rfui+i,«,+2 H 'r du.^i.u, (j - i (mod n) < n/2), 

+ dui_i,ui.2 H ^ duj+i.uj U - i (mod n) > n/2). 



Then, we have 



n—1 n—1 

^ ^ ^ ^u,,u,+k 

vGV v'GV\{v} i^O k^l 

n—1 n—1 

fc=l i=0 



n/2-l„-l 

n—1 n—1 

fc=n/2+l i=0 
n-1 

+ ^ ] (^^U,-,Mi + l + '^tli + l/lli + 2 + • ' • + '^tli + „/2-l/"i + n/2) 
1 = 

^2(l + 2 + ... + (--l)).+ -.— . 

(4) Since u* is a vertex that attains rmiiy^v J2v'£V\{v} dw' , the inequality 
obtained in (3) directly implies the desired one. □ 

Now we discuss the average of the traveling distances of K^^mi''^) ("^ ^ 
{0, 1, . . . , 2ri, — 3}). The traveling distance of a schedule is subject to the follow- 
ing constraint, say the athome constraint: each team stays at its home venue 
before a tournament and returns to its home venue after a tournament. For 
simplicity of the analysis of the approximation ratio, we temporary replace the 
athome constraint with the following assumption. 

Assumption A. If a team plays away games at both the first and last slots, 
then the team moves from the venue of the last opponent to that of the first 
opponent, instead of the moves before the first slot and after the last slot. 

We discuss a traveling distance of each team under Assumption A. Application 
of Assumption A guarantees that a route of each team in i^^RR ('^) ("^ ^ 
{0, 1, . . . , 2n — 3}) is a Hamilton cycle on G (see Figure [2]), and the traveling 
distance of ^^[(^^(to) is invariant with respect to to e {0, 1, . . . , 2n — 3}. Thus, 
we only need to consider K^^^ . This assumption makes the analysis of the 
approximation ratio much easier. 

Let the length of the cycle (wq, wi, . . . , w„_2) obtained by Christofides' 
method in the proposed algorithm be r'. Note that r' < (3/2)r, where r 
denotes the length of a shortest Hamilton cycle on G. Analyzing the structure 
of A' 

DRR reveals the following lemma. 
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Fig. 2 Effect of Assumption A. 



Lemma 4 Under Assumption A, the traveling distance of team t in K^^^ is 
bounded by 



t' + dy^^y. + - '^«t,i-t+i {t e {0, 1, • • ■,n/2 - 1}), 

t' + d„i„i,«* + c?!;',«t - dt;t_i,«t [t e {n/2,n/2 + 1, . . . , n ~ 2}), 

nT/2 (t = n- 1). 



Proof. When t G {0, 1, 2, . . . , n/2 — 1}, team t moves along a Hamilton cycle 
(wt, u*, Wt+i, . . . , u„-2, "^0, wi, W2, • • • , Wt_i). Consequently, the length of the tour 

is T ^ dy^ y* -\- dy* ^y^j^-^ dy^^y^^-^. 

When t 6 {n/2, n/2 + 1, . . . , n — 2}, a tour of team i is a Hamilton cycle 
(wt, ut+i, . . . , u„_2, Wo, wi, W2, ■ • ■ , vt-i,v*), and thus the length is t' + dy^_^,y* + 

Since a tour of team n — 1 is Haniiltonian, Lemma |3I2) implies the desired 
result. □ 



The above lemma implies an upper bound of the traveling distance of Ky)RR- 



Lemma 5 Under Assumption A, the traveling distance of K^)^^ is bounded 
by (71 - 2)t' + 2 + (3/2)r + nT/2. 

Proof. Consider the sum total of upper bounds obtained in Lemma [?] 



{n~ \)t' + L + jit/2 



where 



te{0,l,...,n/2-l} 

+ ^ ^ {dy^_j^,y* +rf^*,Dj —dy^_j^,y^j 

tG{n/2,)i/2+l,...,n-2} 
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It is easy to see that 

\tG{0,l,-:'i/2-l} / \te{l,2,...,n/2} 

^ ^ dvt, vt+i ] + I ^ ^ dyt,v* 

^te{0,l,...,n/2-l} / \te{ri/2-l,n/2,...,n-3} 

\tG{n/2,n/2+l,...,n-2} / \te{ri/2-l,n/2,...,n-3} 

'uey\{u*} te{o,i,...:'i-2} 

< 2 ^ ~ t' + (3/2)t 

i)ev\{i>*} 

where the last mequahty fohows from Lemnia[21[l)- From the above, the lemma 
holds. □ 

Here we drop Assumption A and restore the athome constraint, and con- 
sider the increase of the traveling distance in the following lemma. 

Lemma 6 For each team t, let ^x{t) be the traveling distance oft in K^^^ 
under Assumption A. Then, with the athome constraint the average of the 
traveling distances of team t among K^j^j^(m) (m G {0,1,..., 2n — 3}) is 
bounded by ip^{t) + X]ii'gv\{u} dw' / [n — 1), where v is the home venue oft. 
Proof. For a choice m S {0, 1, . . . , 2ri, — 3}, every team t' different from t plays 
away game with t at first slot just once. Thus, the average length of the moves 
of team t before the first slot is bounded by X]i;'Gy\{i;} dw' / (2n — 2). Similarly, 
the average length of the moves of team t after the last slot is bounded by 
Eu'eyXfu} dvv' /(2f^ — 2). Thus, the average of the traveling distances of team t 
is bounded by £A{t) + Y.v'ev\{v} dvv'/{n - 1). □ 

Summarizing the above lemmas, we have the following theorem. 

Theorem 1 The average of the total traveling distances of schedules K^^^ (to) 
(to g {0, 1, . . . , 2n — 3}) is bounded by 

(71-2)t' + 2 d™. + (3/2)r + nr/2 + ^ ^ dyy,/{n-l). 

v£V\{v*} veVv'GV\{v} 

Lastly we show the approximation ratio of the proposed algorithm. 

Theorem 2 The proposed algorithm is a 2.7 5 -approximation algorithm for 
the UTTP. 

Proof. Let z* be the average of the total traveling distances of schedules 
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iCQj^fj(m) (m G {0, 1, . . . , 2n — 3}). From Theorem [T] and Lemma|21j3)(4), we 
have: 

z*<{n-2)T' + 2 d™. + (3/2)t + 71T/2 + ^ ^ d„,>/{n~l) 

veV\{v} v£V v'£V\{v} 

< (n - 2)(3/2)t + 271T/4 + (3/2)t + ?it/2 + (?iV/4)/(n - 1) 

= (3/2)?ir - 3r + (1/2)7it + (3/2)t + {l/2)nT + {l/A)nT + (l/4)?iT/(n - 1) 

= (ll/4)nr - (3/2)t + (l/4)?ir/(?i - 1) < (ll/4)nr. 

The proposed algorithm output a best of ^drr("^) '= {0, 1, . . . , 2n — 3}), 
and thus the travehng distance of the output is at most z* . Since nr is a lower 
bound of the distance of any double round-robin schedule, this concludes the 
proof. □ 

Let us consider a case that we have a shortest Hamilton cycle H on G. In 
this situation, the following corollary holds. 

Corollary 1 If a shortest Hamilton cycle H on G is given, there exists a 
2. 25 -approximation algorithm for the UTTP. 

Proof. We replace a cycle obtained by Christofides' method in the proposed 
algorithm with a cycle obtained from H by skipping vertex v*. Theorem [T] 
implies that the average of total traveling distances of schedules, say z**, 
obtained by the proposed algorithm is bounded by 

z**<{n-2)T + 2 d,,.,, +{3/2)T + nT/2+Y Y 

v£V\{v*} veV v'£V\{v} 

< nr - 2r + 2nr/4 + (3/2)r + nr/2 + (l/4)nr + (l/4)nr/(n - 1) 
(9/4)nT - t/2 + (l/4)??,T/(n - 1) < (9/4)nr. 

Thus, the approximation ratio is bounded by 2.25 in this case. □ 



4 Computational Results 

In this section, we describe the results of computational experiments using the 
proposed approximation algorithm. 

For the experiments, we took the distance matrices of NL and galaxy in- 
stances from the website [2], because they are the most popular instances 
and one having the largest distance matrix (up to 40 teams) , respectively. We 
ran the proposed algorithm for the UTTP version of these instances; to find 
a short Hamilton cycle, we use Concorde TSP solver [5]. It took less than one 
second to obtain a shortest Hamilton cycle even for the largest case (n = 40). 

To evaluate the quality of obtained solutions, we also tried to find opti- 
mal solutions of UTTP instances with integer programming. Computations 
using integer programming were performed on the following PC: Intel Xeon 
3.33GHz*2, 24GB RAM, Windows 7 64bit, and Gurobi Optimizer 4.5.1 [5] 
with 16 threads as an integer programming solver. For both NL and galaxy 
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Table 1 Results for the UTTP version of NL instances 



n 


approx. 


n * TSP 


gap (%)t 


best UB 
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103978 


23.4 




16 


156828 


119088 


31.7 





tgap is obtained by ( ^'^tsp ~ ^) * ^^^'^ 
* optimal 

*best incumbent solution after 500,000 seconds 



instances: for n = 4,6,8 optimal solutions were obtained; for n = 10, after 
500,000 seconds of computations, branch-and-bound procedures were not ter- 
minated; for n = 12 and larger instances, using integer programming it was 
difficult to find solutions better than those obtained by the proposed algo- 
rithm. 

Tables [T] and [2] show the results of experiments for NL and galaxy instances, 
respectively. The first columns of tables denote the number of teams, n. The 
second ones are the total traveling distance obtained by the proposed algo- 
rithm. The third ones are the value of n times the distance of a shortest 
Hamilton cycle, as a simple lower bound. The fourth ones are the percentages 
of the gap between the second and third columns. 

Like most theoretical approximation algorithms, the obtained gaps are 
much better than the theoretical approximation ratio 2.75 (175% gap). For 
the NL instances and the galaxy instances of up to 20 teams, the gap is 
around 25%. For the galaxy instances of more than 20 teams, the gap is less 
than 20%. Note that the gaps shown in the tables are from the ratio of the 
obtained distance to a lower bound, but not to optimal distance. Therefore 
the gaps between the obtained distance and the optimal value are still better 
than the gaps shown in the tables. 



5 Conclusion 

This paper proposed an approximation algorithm for the unconstrained trav- 
eling tournament problem, which is a variant of the traveling tournament 
problem. The approximation ratio of the proposed algorithm is 2.75, and the 
algorithm yields a solution satisfying the no-repeater and mirrored constraints. 
If a shortest Hamilton cycle on the home venues of the teams is available, the 
approximation ratio is improved to 2.25. Computational experiments showed 
that the algorithm generates solutions of good quality; the gap between the 
obtained solution and a simple lower bound is around 25% for small instances 
(up to 20 teams) and is less than 20% for larger instances. 
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Table 2 Results for the UTTP version of galaxy instances 
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20 
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18.2 




24 
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15.8 
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14.0 
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46410 


12.2 
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55104 
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tgap is obtained by ( ^'^tsp ~ 
t optimal 


1) * 100.0 



*best incumbent solution after 500,000 seconds 
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